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We introduce a hybrid method to determine the neutrino mass hierarchy by simultaneous mea¬ 
surements of responses of at least two detectors to antineutrino and neutrino fluxes from accre¬ 
tion and cooling phases of core-collapse supernovae. The (anti)neutrino-nucleus cross sections for 
®®Fe and ^°®Pb are calculated in the framework of the relativistic nuclear energy density func¬ 
tional and weak interaction Hamiltonian, while the cross sections for inelastic scattering on free 
protons p(t'e,e''')n are obtained using heavy-baryon chiral perturbation theory. The modelling of 
(anti)neutrino fluxes emitted from a protoneutron star in a core-collapse supernova include col¬ 
lective and Mikheyev-Smirnov-Wolfenstein effects inside the exploding star. The particle emission 
rates from the elementary decay modes of the daughter nuclei are calculated for normal and inverted 
neutrino mass hierarchy. It is shown that simultaneous use of (anti)neutrino detectors with different 
target material allows to determine the neutrino mass hierarchy from the ratios of Ve- and t'e-induced 
particle emissions. This hybrid method favors neutrinos from the supernova cooling phase and the 
implementation of detectors with heavier target nuclei (^'’®Pb) for the neutrino sector, while for 
antineutrinos the use of free protons in mineral oil or water is the appropriate choice. 

PACS numbers: 21.10.Gv,21.30.Fe,21.60.Jz,24.30.Cz 


I. INTRODUCTION 

Over the past years considerable progress has been achieved in constraining the mixing parameters in the neutrino 
oscillation framework 0,0, based on various experiments involving atmospheric, solar, and terrestrial neutrinos 0. 
It is now well established that neutrinos have non-vanishing rest masses and that the flavor states Ve, r'/r, and 
are quantum mechanical mixtures of the vacuum mass eigenstates vi, V 2 , and vj, 0. However, currently existing 
data do not determine the neutrino mass hierarchy, l.e., the sign of mass squared difference — Tn\. In 

the case of Am §2 > 0 one refers to the normal mass hierarchy (NH), while Am §2 < 0 corresponds to the inverted 
mass hierarchy (IH). Although a number of techniques has been proposed to resolve the neutrino mass hierarchy, to 
date this question still remains open and presents an important scientific challenge. Recent approaches to resolve the 
neutrino mass hierarchy include methods based on reactor neutrinos 00, various baseline experiment s 0 , Earth 
matter effects on the supernova neutrino signal 0110 , the spectral swapping of supernova neutrino flavors (lll |^ the rise 
time of supernova Pe light curves [10 , the analysis of meteoritic supernova material [l3l| , and detection of atmospheric 
neutrinos in sea water or ice [l0 . 

The neutrino burst from a core collapsing supernova (ccSN) provides valuable Information not only on the explosion 
mechanism, but also about the neutrino mass hierarchy that Is extremely difficult to determine In the laboratory [l0 . 
The first observation of the neutrino burst, from supernova SNI987A in the Large Magellanic Cloud, paved the way 
for neutrino astronomy, and new frontiers in neutrino physics and astrophysics research. [T6l - [l0 . Several supernova- 
neutrino detectors are currently running and a variety of new detectors are proposed to observe galactic ccSN burst [T^ . 
[^ . In the anticipation of the next galactic supernova (the occurrence rate Is Ril-3 events per century fl0). with 
expected high statistics of (anti)neutrino induced detector events, in this paper we introduce a hybrid method to resolve 
the neutrino mass hierarchy, based on neutrino and antineutrino reactions with ®®Fe, ^°®Pb, and with free protons in 
terrestrial neutrino detectors. The aim is to explore novel perspectives for the implementation of Vf, and Pg detectors, 
based on various nuclei and/or protons (in mineral oil and water) as main target material. Since for supernova 
SN1987A mainly the Pe signal has been detected, the role of neutrino Induced detector events for understanding 
the fundamental properties of neutrinos remain largely unknown. While most of the supernova detectors based on 
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free protons or nuclear targets still are primarily sensitive to antineutrinos, the recently developed Helium and Lead 
Observatory (HALO) is sensitive to neutrinos through charged current (CC) interaction mainly with [2lj- 

The hybrid method introduced in this work is based on the simultaneous implementation of two detectors, i.e., one 
for the neutrino and another one for the antineutrino signal from a supernova. Our model calculations show that 
complementary information from two respective detectors can provide constraints for the neutrino mass hierarchy. 

From supernova simulations [22|, [ 2 ^ and the neutrino signal analysis of SN1987A MM we understand that 
(anti)neutrino fluxes detected from a distance have two well distinguished components, i.e., one on a short time-scale 
(^ 1 s) due to emission from the accreting matter onto the proton-neutron star and a long-time scale component 
(Ri 10 s) due to Kelvin-Helmholtz cooling of the nascent proto-neutron star. Thus, being able to trace neutrino and 
antineutrino fluxes from ccSN provide a powerful tool to probe the dynamics of the stellar explosion but also various 
neutrino properties \n\ . In particular, the flavor content of the neutrino signal evolves due to both neutrino 
collective effects and matter effects which can lead to a highly interesting interplay of these and distinctive spectral 
features (Mm. In dense neutrino gas collective oscillations occur due to nonlinear flavor evolution phenomena 
associated to neutrino-neutrino interactions (^ . . The time-dependent fluxes from ccSN can be used to imp rove 

the understanding of the neutrino flavor transformations in regions of high neutrino densities [ 261 - 1 ^ and 

in transitions due to Mikheyev-Smirnov-Wolfenstein (MSW) effects occurring in the matter resonance layers of the 
stellar envelope . More details on the neutrino flavor transformations in supernova are given in recent reviews 

MM; 

A microscopic description of neutrino-nucleus cross sections is crucial for the interpretation of the experimental 
data. It is also important to guide the design of new neutrino detectors with sufficient energy and time resolution 
(and also different thresholds), which are sensitive enough to register charged current events in the energy region 
~ 30 MeV from an intergalactic ccSN. In this particular energy region the largest difference of (anti)neutrino 
spectra between two possible neutrino mass hierarchies is expected, due to the impact of collective and MSW effects 

[iiiliillil- 

Neutrino-nucleus interaction and cross sections in the energy region of supernova neutrinos have already been 
studied within a variety of microscopic approaches, including the shell model |46l - l50l |. the random phase approxima¬ 
tion IRPA I Isil - f^ , continuum RPA (CRPA) [5^4^ , combined CRPA and shell model [5l|, t he Fermi gas 

model (6^65||. auasioarticle RPA [bbl - f^ . projected QRPA and relativistic quasiparticle RPA |72l|. In Ref. 
the cross sections for neutrino-induced particle emissions have been calculated for lead and iron neutrino detectors. 
Event rates from lead-based supernova-neutrino detectors have also been studied in Ref. [73 |. in order to explore the 
prospects for untangling the signatures of various oscillation scenarios. 

In the present study, a microscopic theory framework based on the relativistic nuclear energy density functional 
and weak Hamiltonian is employed in the description of charged-current neutrino-nucleus reactions, including the 
properties of target nuclei, neutrino induced excitations, and weak interaction transition matrix elements [^1^. In 
order to account for the induced events in the detector, the primary particle decay modes of the daughter 

nuclei are described, i.e., 7 decay and emission of one or two nucleons. Neutral current neutrino-nucleus reactions 
have not been taken into account in the present study. Their contributions in the neutron emission channel are small, 
less than r; 10% in heavy nuclei [t^ . To determine positron related signals in water and mineral oil, De cross sections 
of free protons were calculated in the framework introduced in Ref. (ttI . , based on chiral perturbation theory. 

One of the main advantages of the hybrid method established in this work is the simultaneous detection of com¬ 
plementary parts of the neutrino spectra, i.e., Ve and Pe, in at least two different detectors, thus providing additional 
constraints as compared to traditional supernova analyses based mostly on antineutrino events in the detector. The 
time separation and characteristics of the fluxes from different supernova phases are used for the prediction of the 
detector events. Although the absolute number of (anti)neutrino events seen by the detector depends on the distance 
between supernova and Earth, the presented hybrid method is almost independent on this parameter and offers a 
promising approach for solving the problem of the neutrino mass hierarchy. 

The paper is organized as follows. Section |TT] includes the basic formalism and calculations of the (anti)neutrino 
fluxes incoming from an intergalactic ccSN, including the effects that depend on neutrino mass hierarchy. In Sec. 
m we outline the nuclear models used for calculation of exclusive (anti)neutrino-nucleus cross sections and for the 
decay modes of daughter nuclei, i.e., one- and two-particle emissions together with 7 decay. In Sec. IIVI we provide 
further details on the actual hybrid method to determine the neutrino mass hierarchy. The analysis of the detector 
responses to supernova (anti)neutrinos and results of the hybrid method are presented in Sec. |Vl This includes a 
discussion of the role of uncertainties in supernova (anti) neutrino fluxes and their impact on reactions in different 
types of detectors. Section |VT] summarizes the results and includes an outlook for future studies. 
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II. SUPERNOVA NEUTRINOS 


Neutrinos and antineutrinos are produced at various stages of the ccSN evolution [t^ . Model simulations result in 
several characteristic phases of neutrino emission, and provide the evolution of the (anti)neutrino luminosities and of 
the average energies (i) luminosity rise during core collapse, (ii) shock breakout burst, (hi) accretion phase in 

which matter is still falling onto the proto-neutron star (from few tens to few hundreds ms after bounce), ending when 
neutrino heating reverses the infall, and (iv) Kelvin-Helmholtz cooling of the hot proto-neutron star with a duration 
of 10 s or more, accompanied by mass outflow in the neutrino-driven wind. The initial neutrino bursts (i),(ii) are 
unlikely to be used for improved neutrino-signal analysis due to the small number of detector events expected for 
average intergalactic supernova distances (2^. In the focus of the present study are neutrino induced responses from 
(anti)neutrinos associated to the accretion (iii) and cooling (iv) phases of ccSN. 

During the accretion phase, the (anti)neutrino electron species with average energies {E^^) Ri 8 — 12 MeV dom¬ 
inate [IJ m, [13 , with roughly a factor of two larger luminosity than non-electron neutrinos in [13 , while in [13 
the difference is somewhat smaller. Within the cooling phase most of the energy released from the core leads to the 
creation of pairs, and neutrino luminosities may decrease by an order of magnitude. As shown in supernova 
simulation, in the cooling phase the luminosities of electron and non-electron species decrease with time evolution, 
and their differences are small [13 [llj ■ It is important to emphasize that the ordering of average energies remains 
the same through both phases, i.e., {E^^) ^ (Ej^^) ^ {E^^). In this work, however, we will also explore deviations of 
luminosity ratios between different neutrino species to study the respective sensitivity of the detector response in Sec. 

El 


In the inner region of the star, up to several tens of km above the neutrinospheres, neutrino flavor conversion is first 
frozen by synchronization due to the strength of neutrino-neutrino interactions [37, 381. At larger distances from the 


Neutrino-neutrino 


neutrinospheres, the interaction strength decreases and bipolar oscillations occur iimh 

interactions provide a nonlinear term in the equations of motion [13 [11] and the result of such collective effects 
in the high-density neutrino region are spectral splits and swaps in the (anti)neutrino spectra. It is important to 
emphasize that recent studies of neutrino oscillations suggested that collective effects are strongly matter suppressed 
in the accretion phase [sil [slj. In the present analysis we include the three known neutrino species (r'e, and Vt), 
together with their antiparticles. We do not consider sterile types or a fourth generation. 

The matter density profile of the exploding star changes drastically along the neutrino path. At a specific den¬ 
sity, neutrino interactions with matter start dominating through forward elastic neutrino-electron scattering and 
the Mikheyev-Smirnov-Wolfenstein (MSW) effect induces efficient flavor transitions. There are two MSW resonance 
regions, the layer at higher densities (H-resonance layer) and lower densities (L-resonances layer) [4l|. The neu¬ 
trinospheres, however, are far below the resonance layers. In the neutrinosphere, the electron density is orders of 
magnitude lar ger than in the resonance layer. Thus, the flavor transition can simply be described by the Landau- 


Zener formula 142| . In the iron core-collapse supernova the region dominated by collective effects and the two MSW 
resonant regions are spatially well separated (for neutrino energies above few MeV) [13 . Therefore, probabilities for 
flavor transitions in different regions can be factorized [l^j. For the accretion case, the transition probability is given 
by a 3 X 3 matrix, defined as 


p _ pL pH 

^MSW -^MSW 5 


( 1 ) 


due to matter suppression of collective effects, and for the cooling phase we get 

Pc = ^MSW-^MSW-^coU , 


( 2 ) 


where Pcoii represents the flavor transition probability matrix at the end of the collective region, while the matrices 
-^MSW -^MSW related to conversions in the high and low matter density resonant regions, respectively, due to 
the MSW effect. In general, matrix elements are different for neutrinos and antineutrinos, and hierarchy dependent. 
They were calculated in Refs. [13, [U, [13] ■ The arrival of the supernova shock in the outer layers of the star can leave 
a mark on the Ve{i>e) spectra and even cause a non-adiabatic conversion and multiple MSW effects [l^j . In addition 
to a multi-zenith angle instability 1831. i t has recently been shown that a multi-azimuth angle instability can occur 
in the case of normal mass hierarchy |84l| . These effects were not taken into account in the supernova neutrino fluxes 
used here, i.e., a single-angle approximation was employed. 

(Anti)neutrinos leave the surface of a star as mass eigenstates. On arrival at Earth, the incoming neutrino fluxes are 
incoherent due to spread of wave packets for typical supernova distances (~ 10 kpc), thus suppressing oscillations along 
the travelled distance L. In fact, they coincide with the neutrino fluxes on the surface of the star up to a geometrical 
factor 1/A^. Few neutrinos, however, are detected as flavor states. The probability of detecting a particular flavor 
specie is given by 

pEarth _ PpjyffjsPMSW^MSW i (3) 
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for neutrinos emitted during the accretion phase, and 

pEarth _ Pp]y[i,jgP^g^P^g^P(,Qll , (4) 


for emission in the cooling phase. The matrix Ppmns represents absolute squares of the Pontecorvo-Maki-Nakagawa- 
Sakata matrix elements [l|. The above equations are applicable when the detector is directly turned to the neutrino 
source, minimizing Earth effects on the (anti)neutrino spectra. Therefore, the incoming fluxes are given by 




(5) 


where g{r) = is a geometrical factor depending on the supernova-Earth distance r. This is a general expression 

and can even be used with time-dependent (anti)neutrino fluxes. The initial energy distribution for each flavor is 
given by 




( 6 ) 


where is the initial number flux of a particular species, denotes the neutrino 

luminosity, and Ry^ corresponds to the radius of respective neutrinosphere. The initial neutrino spectra can be well 
described using a modified power-law distribution, also known as alpha-fit [s^. 


^^^{Ey^) = N 



{fiy 1) 



(7) 


with the normalization constant 


(/?. + 

(E0)r {Py + 1) ’ 


where T denotes the Euler Gamma function, and fy is the pinching parameter. In general, Py is different for each 
neutrino species and varies in time. This parameter, however, is only weakly constrained by the SN1987A data. Only 
an average value for each phase of spectra can be derived from a fit to the data. Alternatively, theoretical results 
from supernova simulations (including the general time- and phase-dependence behavior) can be used instead. In this 
work, the alpha-fit distribution is employed, rather than the Fermi-Dirac distribution, due to its better description 
of the high energy tails of the spectra 1^. As shown in recent analysis of high-resolution neutrino spectra from a 
spherically symmetric supernova model [87j| . for the purpose of signal forecast in different detectors reasonably good 
accuracy can be achieved by considering the alpha-fit distribution. High energy tails of neutrino distributions are 
important because of the strong increase of neutrino-nucleus cross sections with neutrino energy and the possible 
dependence of neutrino-induced two-neutron emission rates in heavy nuclei on high-energy neutrino spectra. Such 
neutron emissions will be analyzed and discussed in Secs. IIIl| andlVl 
Finally, we also define the relative luminosity as 


ll<^^\At) = 


Lt\At) 


(9) 


in order to avoid an explicit dependence on the total luminosity. In Sec. |V]we explore the sensitivity of the numbers 
of neutrino induced events in detectors and respective quantities introduced within the hybrid method, on the ratio 
of relative luminosities, ly^ : ly^ : ly^. 


III. (ANTI)NEUTRINO-NUCLEUS CROSS SECTIONS 

The charged-current (CC) neutrino- and antineutrino-nucleus reactions are considered in two steps, i.e., primary 
reactions 



(10a) 

h'a + zX N -|- l'^, 

(10b) 
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and decay channels of the daughter nucleus in the case of Eq. (llOal) 

Z+l^N-l + 7 
'z+\^n-2 + n 
2 + P 

A — 4 XA I 

2+i^Af-3 + 2n 

or in the case of Eq. (llObI) 




A 

Z-1 




Xn+1 

-f 

7 

Xn + 

n 


Xn+1 

-h 

P 

Xn-1 

-h 

a 

Xn-1 

-h 

to 


( 11 ) 


( 12 ) 


where la denotes the charged lepton (electron, muon, or its antiparticle). In this paper, however, we are mainly focused 
on electron species. The detector events following charged-current (CC) primary reactions with muon (anti)neutrinos 
are negligible for typical ccSN fluxes. The CC (anti)neutrino-nucleus cross section is given by [s^ 




1 

(^ 


V^PiEi Y, 

lepton 

spins 


1 

2Ji + 1 


E \{f\Hw\i)\\ 

MiMf 


(13) 


where pi and Ei are the momentum and energy of the outgoing lepton, respectively, and E is a normalization volume. 
The Hamiltonian Elyy of the weak interaction is expressed in the standard current-current form, i.e., in terms of the 
nucleon J\(x) and lepton j\{x) currents, as 

Hw = --^j dxJx{x)j^{x) , (14) 

and the transition matrix elements are 

{f\Hw\i) = j dxe~^‘>^{f\J^{x)\i) . (15) 

For the purpose of this work, the cross sections for the charged current zzg(i7g)-nucleus reactions are calculated for the 
following target nuclei: ^^C, ^®Fe, and ^°®Pb. The exclusive cross sections of the primary reactions are described 

in the framework based on the relativistic nuclear e nerg y density functional (RNEDF) [^, [^, [^, by employing 
the density dependent effective interaction DD-ME2 in the particle-hole channel, while pairing correlations are 
described by the pairing part of the finite-range Gogny interaction with set DIS [9^. This framework has been suc¬ 
cessful in the description of Gamow-Teller transitions and other charge-exchange modes, thus it provides a reasonable 
method of choice to calculate relevant neutrino-induced transitions in nuclei 1^. Transition matrix elements 
for neutrino-induced reactions are obtained using the general formalism from Refs. [63l 1^ . This method allows to 
determine the ground state properties of target nuclei and transitions induced by neutrinos in a consistent way. 

The exclusive cross sections are calculated as functions of excitation energy of the initial nuclei, including all 
contributions from the initial ground state of the even-even nucleus to the particular excited state of the daughter 
odd-odd nucleus, for all relevant multipolarities J < 5 and both parities. With the incoming (anti)neutrino fluxes in 
a terrestrial detector, the transition cross section to a particular state i of the daughter nucleus is given by 


{E^JdE, 


^-S'thresh 


(16) 


and the flux averaged cross section is 


(ctz) 


TPtot ’ 

Ua 


(17) 
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where the index i denotes the excited state in the nuclear daughter with excitation energy spin J, and parity 

TT. The quantity i^thresh is the energy threshold for a particular reaction (they usually vary from a few hundred keV 
to « 10 MeV). 

The primary reactions create daughter nuclei which de-excite through emission of photons or particles. Only these 
emissions can be registered in a terrestrial detector. Since we want to explore signatures of neutrino-induced reactions 
in such detectors, the emission channels of the daughter nuclei, as shown in Eqs. m and () 12 | 1 . have to be followed in 
a calculation. In this calculation we included all the channels shown in the above equations but here we particularly 
focus on neutron emission, as charged particle emission is negligible in comparison. The neutrino flux-integrated 
emission cross section ai^ = criPi^ for daughter state i and a particular channel k is simply the product of the 
flux-integrated cross section ai and a branching ratio Pi^, describing the relative contribution of channel k within all 
possible decay channels. This branching was calculated from the ratio of transmission coefficients obtained from 
solutions of the time-independent Schrddinger equation in the optical model, 

P,= ,K,K'^2n. (18) 


These transmission coefficients are related to reaction widths of resonances and are similar to those used in the 
Hauser-Feshbach model which describes compound nucleus reactions Therefore, we used the Hauser-Feshbach 
code SMARAGD 941 to determine Pi^- The procedure to calculate transmission coefficients is similar to the 
one described in |94l . l95l| . Emissions to final states include emissions to experimentally known discrete levels up to a 
final state excitation energy Egn. Above Egn an integration over a nuclear level density for spins J < 12 and both 
parities is invoked. Thus, the calculated emission spectrum contains both discrete and continuous contributions, as 
is also typical for compound nucleus reactions. Here, the “compound” nucleus is formed by the primary CC weak 
reactions. The same default settings of the code were used as in [9^ but with updated experimental levels (up to 
40 known excited states from the 2010 version of NuDAT [9^) and level density (including parity dependence as 
described in [^1. Emission thresholds are given by the particle separation energies in the emitting nucleus. These 
were computed from measured nuclear mass differences taken from . It should be noted that the 7 cascade in the 
daughter was not followed in detail and therefore neutron emission after 7 deexcitation(s) is not included. We expect 
it to be negligible, anyway, because significant contributions to the total particle emissions only come from at with 
moderate excitation energy, < 25 MeV, and separation energies are of the order of 10 MeV. 

For emission of the second neutron, Tin 2 = was calculated from a recursive procedure, starting 

from the transmission coefficient for one-neutron emission from state i to final state j. Tin. For all final states j 
possibly populated by emission of one neutron, the procedure of determining a branching ratio P* for single-neutron 
emission from the new daughter nucleus in state j was repeated. 

Using the above approach, it is not only possible to calculate flux-integrated cross sections for the individual particle 
emission channels but it is also possible to determine the emission spectrum ^^(Fem), he., the number of particles 
emitted with a given energy Fem- This is useful if the detector cannot only provide an energy-integrated signal but also 
resolve the actual particle spectrum. To illustrate a first application of the hybrid method, however, in the following 
we restrict ourselves to predicting particle counts in the detector, i.e., energy integrated results. The detailed emission 
spectra are available from the authors on request and may provide further constraints in future analyses of neutrino 
signatures in detectors. Since the studies of lead detector efficiency are still in progress in the present analysis 
perfect detector response is assumed, with the efficiency denoted by ry —> 1 . 

The rate of particles registered in the detector is given by 


dN^ 

dt 


= VtUk ^ CTiPiK , 


(19) 


where in the case of emission of only one neutron (In) from daughter nucleus we have Pun = Pun — Pam while for 
7 , one-proton (Ip), two-neutron (2n) and emission of a particle = Pin- The total number of targets (nuclei or 
protons) in the detector is = Xm^Nj^/M,., where m<i denotes the detector mass, X the number of particular 
targets (e.g., X = 2 corresponds to the number of protons in water molecule, otherwise if no molecules are formed, 
X = 1), Mr the respective molar weight, and Na the Avogadro constant (see Sec. IIVI for more details). 

In addition to z/e(i^e)-nucleus cross sections, in this work also inelastic 9e scattering on free protons is taken into 
account. This is relevant for detectors using mineral oil (alkanes, which contain CH 2 groups) or water (H 2 O). In the 
calculation, we employ heavy-baryon chiral perturbation theory, which also includes radiative corrections lz3,[2ii- 
The calculated cross sections for inelastic De+P scattering are shown in Fig.[TJ In Sec El this cross section is employed 
in modeling the number of the detector events induced by antineutrinos in mineral oil and water. 
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IV. DESCRIPTION OF THE HYBRID METHOD 


Predictions of the absolute values of the number of events for NH or IH may not provide conclusive information on 
the neutrino mass hierarchy due to the overlaps of the expected value ranges caused by the considerable uncertainties 
in the modeling of supernova (anti)neutrino fluxes [^ . In the following we introduce a hybrid method that could 
constrain the neutrino mass hierarchy. It is based on observing CC reactions with at least two detectors which cover 
complementary parts of supernova (anti)neutrino spectra, i.e., one detector is sensitive mainly to Vg while the other 
to Pe- At the core of the hybrid method is the definition of a set of hierarchy-dependent quantities, involving ratios 
of event numbers from several detectors and making use of different reactions as well as different phases of supernova 
evolution. In this way the dependence on the absolute values of some of the involved quantities is reduced, e.g., on the 
luminosities, on the supernova-Earth distance, and on the mass of the detector. In this work, we explore the feasibility 
to provide conclusive information on neutrino mass hierarchy by using the relative quantities of such a hybrid method 
based on two detectors, rather than depending on the absolute number of the events in a single detector. 

In the analysis of time-dependent detector events, for each reaction channel the ratio of the numbers of emitted 
particles, both for the supernova accretion and cooling phase, is given by. 


TD&{c) ;AB 
^kai^b 


(At,) 




( 20 ) 


and similarly its time-integrated variant. The indices A and B identify the detector with its particular decay or 
deexcitation channel ka, kb, of the daughter nucleus (see Sec. Ell. In addition, the labels ta and tb identify whether 
nuclei or free protons are involved in CC reactions in the detectors. For example, there are several different target 
nuclei in a water or mineral oil detector (free protons and isotopes of C or O), which have different contributions to 
the number of CC events due to their different CC reaction cross sections. We assume that the signal of a particular 
emission mode can be well distinguished from signals of other modes. Furthermore, the time-bin Atj should be 
wide enough to allow for good statistics of the events and it should be similar (comparable) for both detectors. 
If the v-v interaction is sufficiently strong to cause collective effects in the region up to few hundred km above 
neutrinosphere, then flavor transitions will occur in complementary parts of the spectra, causing spectral splits and 
swaps. Furthermore, in the case of NH the MSW effect in both resonant regions will contribute to transitions between 
neutrinos, while the case of IH is characterized by transitions in the lower resonant region only. This is different, 
however, for antineutrinos which experience transitions in the low resonant region for NH and in both resonant regions 
for IH (for details, see [llj). The quantity is a relative measure of total supernova effects on initial (anti)neutrino 
spectra. 

From Eqs. m and dsni) we get 


<l^''(At,)=AAB 


TA 

iKA 


D 




'TB 

IKB 


where 


Aab 


rriANA 

Mfi- 

tubNa ’ 


( 21 ) 


( 22 ) 


is constant for a given combination of neutrino (B) and antineutrino (A) detector. The coefficient stands for 

the number of targets (atoms or protons) in a molecule for mineral oil or water detectors. In the case of a detector 
using a metal, = 1. The relative molecular (A) or atomic (B) molar weight is denoted by Note that 

even in a single detector several targets can give contributions to the total 7 , In, Ip and 2n emissions. In Eq. (ED 
one can see that R^^k,b contains three types of information, the coefficient Aab depends only on the combination of 

detectors used, the coefficients Dk and are connected to the particle emission and its probability of observation, 
and third quantity aj, which is the transition probability at a given neutrino flux per target nucleus. The latter term 
is rather complex as it depends on the structure of the nuclear ground and excited states, and on the (anti)neutrino 
spectra. Thus, it also contains imprints of collective and MSW effects, which are generally different for NH and IH. 

In order to eliminate the dependence on the particular detector combination, thus using only information on 
particular targets involved in CC reactions for 77 —>■ 1 (perfect detector response), we define the reduced quantity, 


= lim R 


AB 

KaKb 


/Aab , 


T 


(23) 







which is an effective ratio of transformation probabilities (including decay modes) per single atom (B) or molecule 
(A) caused by complementary parts of the incoming supernova neutrino fluxes. The quantity still contains 

information about targets included in the processes, and on the empirical structure for molecules when dealing with 
water and mineral oil. In other words, the reduced quantity is good measure for hierarchy-dependent effects, 

i.e., a relative measure of the effects that shape (anti)neutrino spectra, independent of the absolute values of detector 
mass, effective number of targets, and explicit supernova-Earth distance. 


V. RESULTS 

The results of the present analysis are twofold: the first category is related to the detector response to the fluxes of 
supernova (anti)neutrinos, i.e., the induced particle emissions from target nuclei; the second is related to applications 
of the hybrid method to constrain the neutrino mass hierarchy. 


A. Supernova (anti) neutrino fluxes 

In order to investigate the feasibility of the hybrid method despite of the current uncertainties in predicting neutrino 
fluxes [13, the effects of various sets of initial supernova flux parameters are explored. We used several combinations 
of initial luminosities in the accretion 10®^ erg/s) and cooling phase 10®^ erg/s) fluxes, with a fixed total 
emitted energy for each phase (see Tab.|T|, i.e., Ri 0.4 x 10^^ erg for a duration of 0.4 s, and Ll°°^ k.2x 10®^ erg 
for a duration of 10 s. We also varied the initial average energies of neutrinos from 8 to 12 MeV for from 11 to 
15 MeV for and from 11 to 19 MeV for non-electron species. Thus, the range of values from different simulations 
of supernova neutrino fluxes is covered (lH . [lH , and the sensitivity of the hybrid method on the differences in 
neutrino fluxes can be explored. In the present analysis, the initial averaged energies of neutrino species were kept in 
the canonical order, as shown in Sec. nil with two additional conditions: {EpJ- < 3 MeV and {E°J- (E,° ) <4 
MeV. The pinching parameter in Eq. © was set to 4.0 for the accretion and 3.0 for the cooling phase fluxes [^.1^. 
In calculating transition probabilities beyond the collective region, we have used fixed values of splitting energies from 
a supernova neutrino flux analysis [s^ - The incoming (anti)neutrino fluxes were calculated including collective and 
MSW effects in the core-collapsing star (see Sec. II). Model calculations include best-fit values o f neutrino oscillation 
parameters. They can be found, e.g., in [43, although somewhat different values are given in [lOOl . Il0lj |. As a test 
case for ccSN in our galaxy, an imaginary star 25000 l.y. away from Earth was assumed. 

As an example. Fig. shows the set of incoming Vg and 9^ fluxes of the accretion and cooling supernova phases 
for the normal (NH) and inverted mass hierarchy (IH). Three combinations of the luminosity ratios of (anti)neutrino 
species for type II supernova spectra in accretion(Al-A3) and cooling (C1-C3) phase are considered, as given in Tab.lH 
The initial average energies of the supernova fluxes shown in Fig. [5] are {E^ ) = 10 MeV, {E^ ) = 13 MeV, and {E^ ) 
= 15 MeV. 


B. Nuclear responses to supernova (anti)neutrinos 

Figures [3] and |4] show examples of the CC primary cross sections of ®®Fe averaged over supernova cooling phase 
fluxes of t'e and 9e (for the case C3 in Tab. |T|, respectively, and displayed as a function of excitation energy in the 
initial nucleus. Contributions to the cross sections are shown separately for several multipoles, from J'^ = 0^ to 
=3^. The energy threshold for CC reactions in ^®Fe is relatively low for Ve (4.57 MeV) and Ve (4.72 MeV). The 
z/e-^^Fe cross section is mainly determined by = 0+, I’*' transitions, while the contribution of = 1~,2~ states 
is an order of magnitude smaller. Other multipolarities contribute only marginally. As shown in the upper panel 
of Fig. m the z7e-^®Fe reaction cross section, on the other hand, is dominated by = I’*' transitions, with a major 
contribution from the excited state at 4.5 MeV. The lower panels of Figs. |3] and |4] show the ratio of flux averaged cross 
sections between two hierarchies for the specific incoming cooling phase Vg or Vg flux C3. They indicate that exclusive 
transitions to particular states should be reduced in the case of IH by factors ~ 1.1 — 2.15 for Ue and ~ 1.15 — 1.32 
for 9e. A similar effect is expected for the other flux cases. 

The total z/e-nucleus CC reaction cross section is increasing rapidly with increasing number of neutrons in the 
nucleus [73. The relative difference in the depths of the proton and the neutron part of the mean field potential 
increases in nuclei with neutron excess, however, Pauli blocking strongly suppresses transitions from proton to neutron 
quasiparticle states in Pe-nucleus reactions. This effect is evident in the Pe-^^Fe reaction, where the cross sections are 
an order of magnitude smaller than in the case of neutrino induced reaction (see Figs. [3]and|4|), although ®®Fe has 
a relatively small neutron excess {N — Z = 4). The upper panel of Fig. [5] shows the exclusive cross sections for the 
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j^^_208pb reaction, illustrating an example of a very neutron-rich target. The overall cross sections are considerably 
larger than for ^®Fe, thus indicating that a lead-based detector is a reasonable choice to efficiently detect supernova 
neutrinos. On the other hand, due to the strong blocking effect in the single-particle spectra, i/e-^°®Pb cross sections 
are considerably reduced, i.e., the lead detector allows measurements of supernova-neutrino induced events only. 


C. Detector responses to supernova (anti) neutrinos 

By employing the framework outlined in Sec. nil we study neutrino induced events in detectors based on various 
target material. Since perfect efficiency is assumed, we note that in the realistic case the observed number of the events 
may considerably redu ce, e. g., the efficiency to detect a neutron generated in a lead-based detector can be as low as a 
few tens of percent Il02l |. Furthermore, it is assumed that the detector is turned directly to the incoming neutrino 
flux, thus eliminating Earth effects on neutrino spectra. Four cases of target material were considered, mineral oil 
(CH 2 ), water (FI 2 O), ®®Fe, and ^°®Pb. The mass of the detector was taken to be 1 kt for all cases of target material. 

The numbers of (anti)neutrino-induced events due to In, 2n, and emissions for several different accretion 
supernova fluxes are shown in Tab. m both for normal and inverted neutrino mass hierarchies. Three combinations 
of the luminosity ratios (A1,A2,A3) are considered, as given in Tab. H] In addition, three sets of the initial average 
energies of the supernova fluxes are used, (i) {E^J = 8 MeV, {E^J = 11 MeV, {E^J = 13 MeV; (ii) {E^J = 10 
MeV, {E^J = 13 MeV, {E°J = 15 MeV; and (hi) = 12 MeV,° = 15 MeV, = 19 MeV. The number 
of i/e-related events, i.e., the number of emitted positrons, is expected to depend not only on the initial average 
energies of the z/ spectra, but also on the initial luminosity ratios between neutrino species. Due to large overlap of 
the respective positron events, however, it is difficult to put any constraint on the neutrino mass hierarchy from these 
events alone. 

The number of single-neutron emissions in ®®Fe (^°®Pb) is ~ 10 (few) times smaller than positron events related 
with free protons in water and mineral oil, while two-neutron events in ^°®Pb are negligible for cases (i) and (ii). The 
number of single-neutron events in ®®Fe and ^'^^Pb significantly increases with the initial average energies of the ly 
spectra. The number of emitted neutrons, however, remains ambiguous for the shown cases, due to the strong matter 
suppression of collective effects, which was also observed in [s^]. Therefore, the analysis of the detector response to 
the (anti)neutrinos from the accretion phase does not allow to distinguish between the two neutrino mass hierarchies. 

Table nni shows the number of detector events in the same target materials as discussed above, but for the 
(anti)neutrino fluxes from the cooling phase, for the luminosity ratios in Tab. U denoted as C1,C2,C3. Again, three 
sets of the initial average energies of the supernova fluxes are considered. The number of positron events in water 
or mineral oil induced by is at least a few times larger than the number of primary Vg reactions with ®®Fe and 
208pb fQj. all investigated cooling phase neutrino fluxes (with different values of initial average energies as shown in 
the table). The number of single-neutron emissions for ®®Fe and for one- and two-neutron emissions for ^*^®Pb target 
is larger in the case of NH (see Figs. [3] and [S]). Contrary to the accretion phase, the neutrino fluxes from the cooling 
phase result in an increased number of the events for NH and thus in a reasonable separation between the numbers 
of the particle emissions from nuclei for the two hierarchies. 

The most pronounced response to CC interactions with is obtained for ^°^P b. Due to large number of protons 
{Z = 82), Coulomb effects enhance the phase space for emitted electrons |l03j |. It also has large neutron excess 
{N — Z = 44) and relatively low threshold (2.88 MeV) for CC reactions. For the supernova cooling phase, we obtained 
a flux averaged cross section for neutron emission from of ~ 10“"^° cm^. In addition, a significant contribution of 
two-neutron events was obtained due to relatively high two-neutron cross sections (~ 10 “'^^ — 10 “^° cm^) at energies 
above the two-neutron separation energy of 23 MeV relative to the mother nucleus (see upper panel of Fig. [5|). The 
same applies for all sets (i)-(iii) of neutrino fluxes. The sensitivity of the number of 2n-emissions on the type of mass 
hierarchy is evident in the cooling phase, i.e., the ratio Ri 2 (see lower panel of Fig. [5] and Tab. IIIII) . 

The cross sections for i^e—^°®Pb reactions are strongly suppressed (by three orders of magnitude) due to Pauli 
blocking of the neutron single-particle states. Thus, lead presents an excellent choice for a supernova neutrino 
detector to be used in the hybrid method. In order to obtain good statistics of the events, though, a sufficient amount 
of target material has to be available, e.g., Ri 1 kt when assuming typical intragalactic supernova-Earth distances. 
The few reactions with be and in this material can be neglected in comparison to related events. Moreover, 
natural lead has a comparatively small thermal neutron capture cross section of r; 0.1 5 ba rn, which allows lead to 
moderate fast neutrons almost without absorption prior to arrival at a neutron counter |l04| . Althou gh le ad isotopes 
have large neutron excess, some of the neutrons might still be captured by ^°®Pb and ^°^Pb, though |l05j |. 

In comparison, a ®®Fe target is characterized by somewhat weaker response to CC reactions with neutrinos (see 
Tab.lU and Tab. nni for the accretion and cooling phase neutrinos, respectively). The predicted total cross section 
of neutrino-induced reactions in ®®Fe is r; 10 times smaller than for ^°®Pb. Almost half of the predicted events come 
from single-neutron emission from ^®Co, while two-neutron emissions are severely reduced. Due to small neutron 
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excess {N — Z = 4), the detector based on ®®Fe is also sensitive to antineutrino CC reactions with cross sections of 
Ri 10-42cni2^ order of magnitude smaller than for z/g—®®Fe reactions, but the daughter nuclei mostly deexcite by 
emission of 7 rays. Furthermore, the cross sect ions for absorption of thermal neutrons in natural iron is an order of 
magnitude larger than in lead (r; 2.5 barn [l04| l. and thus iron appears to be an inappropriate target material for a 
t'e and 9e detector. 

For mineral oil (water) with a density of 0.85 g/cm^ (1.0 g/cm^ at 4 °C), in which the target nuclei are (^®0), 
the predicted cross sections are of the order of R 10“"^^ cm^. Due to iV = Z = 6 {N = Z = 8), the difference in 
total number of primary events between neutrino and antineutrino CC reactions is smaller than for N > Z nuclei. 
Due to relatively high energy thresholds (^ 11 MeV) for both types of reactions in ^^C (^®0), the expected responses 
are rather low for (i/g) energies of < 20 MeV. Actually, in these detectors Pg induce reactions mainly with free 
protons (see Tab. m and Tab. EH, with a low energy threshold (1.8 MeV) and cross sections ct R 10 '^^(Aj^/MeV)^. 
The dominance of z^g-p reactions in mineral oil and water ensures an efficient coverage of the z/g spectra. Therefore, 
detectors based either on mineral oil or water can cover the antineutrino response in the hybrid method. In these 
cases the influence of CC reactions of I'e and Pe with nuclei amounts only r 1% in the response and thus CC events 
with nuclei can be neglected. Similar results were obtained in Ref. (lOOj . i.e., only r 1.5%(5.6%) of the reactions 
are expected to be related to ^®0(^^C) in CC and NC reactions in the KII detector (Large Volume Detector) in the 
energy range of interest. An additional 7 signal related to neutron capture in the detector (either on the original 
target nuclei or after adding a n im purity such as Gd to increase the capture rates) can be used to confirm or eliminate 
charge exchange i^g-p events (l07l| . 


D. Application of the hybrid method 

In the following we present calculations of the ratios of Pe and Ve detector events for a variety of possible incoming 
supernova (anti)neutrino accretion and cooling phase fluxes (see Tab. H]), both for normal (NH) and inverted (IH) 
mass hierarchy. By using as defined in Eq. (1^51) . we obtained mass and distance invariant values for the most 

sensitive decay channels in the detectors. Statistical uncertainties in , however, still depend on such quantities 

as distance, etc. Due to the reasons given above, we eliminated ^®Fe from the further analysis and used only ^°®Pb 
and mineral oil (water) for the neutrino and antineutrino sectors, respectively. Including only dominant channels in 
the analysis, the quantity should be the same in the Pb/H 20 and Pb/CH 2 detector combinations, due to the 

similarity of structural (empirical) formula of water and mineral oil with 2 free protons each. 

In Fig. [5] the following ratios of the detector events are shown for the incoming (anti)neutrino fluxes of 

the accretion phase for NH and IH, including statistical uncertainties: (a) (b) and (c) 

The quantities (a), (b), and (c) correspond to the ratios between the numbers of positron events in the antineutrino 
detector (with dominant Pg-p events) and neutrino-induced one-, two-, and total neutron emissions in the Pb detector, 
respectively. We employed three combinations of the luminosities for the accretion phase (A1-A3) for three different 
configurations of initial average energies of neutrino spectra, as specified in Tab. HI 

As seen in Fig. the ratios of r'g and z/g induced events result in reasonable separation between NH and IH for the A2 
case, which has a larger asymmetry in relative luminosities between electron and non-electron species, independent of 
initial average neutrino energies. The other cases of neutrino fluxes are characterized by partial overlaps between NH 
and IH event ratios for all energy configurations. By considering other ratios shown in Fig. | 6 l we note that neutrinos 
from the accretion phase systematically do not provide conclusive information to determine the mass hierarchy, with 
the exception of flux A 2 . 

In Fig. [7] we show the same decay channel analysis as in Fig.jBJ but for the neutrino fluxes of the supernova cooling 
phase. Three types of the cooling phase (anti)neutrino fluxes were considered (C1-C3), as given in Tab. HI and the 
same energy configurations as for the accretion fluxes given above. As seen in the figure, the (anti)neutrino induced 
events for the cooling phase are characterized by complete separation of the ranges of values calculated for NH and 
IH: (a) (b) (c) 4 '^™®^^for all cases of supernova neutrino fluxes. Since the detectors with 

target material based on heavy nuclei are usually used only as counters, the number of emitted neutrons can provide 
additional information for the zzg part of spectra. Nevertheless, the strong dependence of e"*", one- and two-neutron 
emissions on the average initial energies of v can be used to constrain their initial energy configurations (see Tab. HII|. 

The results of the present analysis also display sensitivity to the degree of asymmetry of the initial relative lumi¬ 
nosities between zz species. In the case of C2 (larger G asymmetry), the r-quantities are systematically shifted toward 
smaller values (see Fig. [7|). In the case of two-neutron emissions, a complete distinction of the quantities 

between two hierarchies is obtained, similar as in the case of one-neutron emissions. Therefore, in a detector with 
sensitivity to one- and two-neutron emissions, such as HALO, both types of the events can be used to constrain the 
neutrino mass hierarchy. We note that while this may be feasible for the (ii) and (iii) cases of supernova neutrino 
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energy configurations, in the (i) case the number of two-neutron events is rather small (see Tab. III). Thus, it would 
be difficult to provide reliable limits. Nevertheless, the separation of the r-quantities as shown in Fig. [71 especially in 
the case of a reasonable number of single-neutron and total neutron emission events (Tab. Ill), clearly demonstrates 
the promise in using the hybrid method to constrain the neutrino mass hierarchy. Additional, currently not available, 
information on the sensitivity and efficiency of the HALO detector would allow a further, more detailed feasibility 
study for the proposed method. 


VI. CONCLUSION 

We have presented a hybrid method to determine the neutrino mass hierarchy by simultaneous measurements of 
supernova Ve and 9e events in detectors based on different types of target material. Using supernova r'e(i^e) fluxes 
that include collective and MSW effects for the accretion and cooling phases, responses in mineral oil, water, ^®Fe and 
208pb iiave been analyzed both for normal and inverted neutrino mass hierarchies. The analysis of charge-exchange 
excitation spectra for ®®Fe and ^°®Pb targets demonstrates the sensitivity of the nuclear response to the neutrino mass 
hierarchy, that is different for various multipole transitions. The hybrid method, that combines antineutrino-induced 
events in water or mineral oil with neutrino-induced emissions in heavier target nuclei (such as ^°®Pb), could provide a 
useful tool to constrain the neutrino mass hierarchy. The number of emitted neutrons for NH in general is larger than 
for IH, both for ^°®Pb and ^®Fe target nuclei. Since it is rather difficult to compare the absolute values of calculated 
particle emissions with the detector events, a set of hierarchy dependent relative quantities is introduced, which are 
independent of the supernova-Earth distance, of the effective number of targets in r'e and Pg detectors, of their total 
masses, and of the ratios of these quantities. They contain the information on nuclear structure, excited states, 
relevant charge-exchange transitions and particle emission modes induced in nuclei by supernova (anti)neutrinos. 
The incoming (anti) neutrino fluxes employed in this study include the imprints of collective and MSW effects during 
the ccSN event and these are generally different between NH and IH. Model calculations show that the ratio of the 
numbers of Vg and Pg induced detector events represents a quantity that allows to distinguish between the two neutrino 
mass hierarchies. Simultaneous comparison of several ratios of the numbers of detector events associated to various 
emitted particles for different supernova phases would reduce the uncertainties that may arise by considering only a 
single quantity. 

As confirmed in calculations, heavy nuclei are almost completely inert to CC reactions with antineutrinos due 
to the strong Pauli blocking of neutron single-particle states. The pure neutrino signal and large CC cross section 
for reactions with neutrinos establish ^*^®Pb as probably the most important nuclear target for neutrino detection 
and reconstruction of the neutrino part of spectra. On the other side, detectors based on mineral oil or water as 
target-material represent the most reasonable choice for Pg detection. Takin g int o account current neutrino detector 
developments (in particular, HALO (^°®Pb) [2l|, Super-Kamiokande (water) |l08| . etc.), the hybrid method presented 
in this work will provide a useful tool to constrain the neutrino mass hierarchy when the next galactic ccSN appears. 
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TABLE I: Initial lumininosity ratios of (anti)neutrino species for type II supernova spectra in accretion(A) and cooling(C) 
phase. Three possible combinations of the luminosity ratios for (anti)neutrino species are shown for each phase, together with 
fixed total energy emitted by (anti)neutrinos and time duration of the phases. 




Tint 

Accretion ph. 

« 0.4 

X 10®® erg 


A1 

4.0 : 

4.0 

: 3.0 

0.4 s 

A2 

3.0 : 

3.0 

: 1.0 


A3 

12.0 : 

8.0 

: 5.0 


Cooling ph. 

« 2 X 

10®® 

erg 


Cl 

1.0 : 

1.0 

: 1.0 

10.0 s 

C2 

2.0 : 

2.0 

: 4.0 


r’Q 

Q O . 

O Q 

. Q c: 




E- [MeV] 

FIG. 1: The cross sections for the inelastic scattering of on free protons, shown as a function of the incomin g an tineutrino 
energy. Calculations are based on heavy-baryon chiral perturbation theory with radiative corrections from Ref. [7^ . . 
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TABLE II: Detector response for i/e(t'e)-induced reactions in mineral oil (CH 2 ), water (H 2 O), ®®Fe and ^°®Pb, for the incoming 
(anti)neutrino fluxes of the accretion supernova phase, both for normal (NH) and inverted (IH) neutrino mass hierarchies. Only 
dominant emission channels are shown, including e^, one- and two-neutron emissions. Al, A2, A3 denote the combinations of 
the luminosity ratios of (anti)neutrino species given in Tab. [H while (i) stands for (IfSe) = 8 MeV, {E°^) = 11 MeV, {E^^) = 
13 MeV; (ii) is = 10 MeV, (Sgj = 13 MeV, = 15 MeV; and (iii) is {El^) = 12 MeV, (ElgJ = 15 MeV, = 

19 MeV. 


Accretion 

ph. 

L type 


(i 

i) 



(ii) 



(iii) 


0.4 s 



NH 


IH 

NH 


IH 

NH 


IH 

p(i5e,e+)n 

in CH 2 














e+ 


Al 

45 

- 60 

42 

- 56 

54 

- 70 

49 

- 64 

65 

- 82 

61 

- 78 



A2 

58 

- 74 

28 

- 39 

69 

- 87 

32 

- 45 

82 

- 101 

41 

- 55 



A3 

43 

- 57 

34 

- 47 

51 

- 67 

40 

- 54 

61 

- 78 

51 

- 66 

p(i5e,e+)n 

in H 2 O 














e+ 


Al 

34 

- 47 

32 

- 44 

41 

- 55 

37 

- 50 

49 

- 65 

47 

- 62 



A2 

44 

- 58 

21 

- 31 

53 

- 69 

25 

- 36 

62 

- 79 

31 

- 43 



A3 

33 

- 45 

26 

- 37 

39 

- 53 

30 

- 43 

47 

- 61 

39 

- 52 

56Fe(:/e,e- 

■)56Co 














In 


Al 

0 - 

4 

0 - 

3 

2 - 

6 

2 - 

6 

6 - 

13 

6 - 

13 



A2 

0 - 

2 

0 - 

2 

1 - 

4 

1 - 

4 

4 - 

9 

4 - 

9 



A3 

0 - 

3 

0 - 

3 

1 - 

5 

1 - 

5 

5 - 

11 

5 - 

11 

In 


Al 

12 

- 20 

12 

- 20 

19 

- 29 

19 

- 29 

35 

- 48 

35 

- 48 



A2 

7 - 

14 

7- 

14 

12 

- 20 

12 

- 20 

22 

- 33 

22 

- 33 



A3 

10 

- 17 

10 

- 17 

16 

- 25 

16 

- 25 

29 

- 41 

29 

- 41 

2n 


Al 

0 - 

3 

0 - 

3 

1 - 

5 

1 - 

5 

7 - 

14 

7 - 

14 



A2 

0 - 

2 

0 - 

2 

0 - 

4 

0 - 

4 

4 - 

10 

4 - 

10 



A3 

0 - 

2 

0 - 

2 

1 - 

5 

1 - 

5 

6 - 

12 

6 - 

12 

total n 


Al 

15 

- 24 

15 

- 24 

26 

- 38 

26 

- 38 

56 

- 72 

56 

- 72 



A2 

9 - 

17 

9 - 

17 

17 

- 26 

17 

- 26 

36 

- 49 

36 

- 49 



A3 

12 

- 20 

12 

- 20 

22 

- 32 

22 

- 32 

46 

- 61 

46 

- 61 
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TABLE III: The same as Tab. El but for the cooling supernova phase. Cl, C2, C3 denote the combinations of the luminosity 
ratios for (anti)neutrino species given in Tab. [J 


Cooling ph. 


L type 


(i 

0 



(ii) 



(iii) 


10 s 




NH 

IH 

NH 

IH 

NH 

IH 

p(t'e,e+)n 

in 

CHa 














e+ 



Cl 

253 

- 285 

248 

- 281 

303 

- 339 

295 

- 331 

394 

- 435 

370 

- 409 




C2 

217 

- 248 

211 

- 242 

270 

- 304 

257 

- 291 

378 

- 418 

340 

- 378 




C3 

235 

- 267 

230 

- 262 

286 

- 320 

276 

- 310 

381 

- 421 

352 

- 391 

p(t'e,e+)n 

in 

H 2 O 














e+ 



Cl 

195 

- 224 

192 

- 220 

234 

- 266 

228 

- 259 

305 

- 341 

286 

- 321 




C2 

167 

- 194 

163 

- 190 

209 

- 238 

199 

- 228 

292 

- 327 

263 

- 296 




C3 

181 

- 209 

178 

- 205 

220 

- 251 

213 

- 243 

294 

- 330 

272 

- 306 

®®Fe(i^e,e- 

)®«Co 














In 



Cl 

11 - 

19 

7 - 

14 

21 - 

31 

13 - 

21 

46 - 

60 

27 - 

39 




C2 

14 - 

23 

9 - 

16 

25 - 

36 

15 - 

24 

55 - 

71 

30 - 

42 




C3 

12 - 

20 

7 - 

14 

22 - 

32 

14 - 

22 

48 - 

63 

28 - 

40 


2°®Pb(i/e.e- 

-)208Bi 







In 

Cl 

88 - 108 

68 - 85 

126 - 150 95 - 116 

194 - 223 

139 - 163 


C2 

106 - 128 

82 - 101 

151 - 177 

111 - 133 

230 - 262 

154 - 180 


C3 

93 - 113 

71 - 89 

133 - 157 99 - 121 

204 - 234 

144 - 169 

2n 

Cl 

9 - 16 

4 - 9 

19 - 29 

10 - 17 

53 - 69 

28 - 40 


C2 

11 - 18 

5 - 11 

23 - 34 

11 - 19 

63 - 80 

31 - 44 


C3 

9 - 16 

4 - 10 

20 - 30 

10 - 18 

56 - 72 

29 - 41 

total n 

Cl 

112 - 134 82 - 101 

173 - 201 

122 - 145 

313 - 350 

205 - 235 


C2 

135 - 159 

98 - 119 

207 - 237 

141 - 166 

371 - 410 

228 - 259 


C3 

118 - 141 

86 - 105 

182 - 210 

128 - 151 

329 - 366 

213 - 243 
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E [MeV] E [MeV] 


FIG. 2: Incoming (left side) and Ve fluxes (right side) for the accretion (A) and cooling (C) phase of core-collapse supernovae 
as a function of neutrino energy for normal (NH) and inverted neutrino mass hierarchy (IH), for initial v energy configuration 
(ii). Three combinations of the luminosity ratios for (anti)neutrino species used in the accretion (Al-3) and cooling (Cl-3) 
phases are given in Tab. [H 
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FIG. 3: Multipole decomposition of the flux averaged cross sections of the i/e-®®Fe reaction as a function of excitation energy 
in the initial nucleus, calculated for the incoming cooling phase neutrino spectra (C3) in normal hierarchy (NH), including 
multipoles = 0^ — 3^ (upper panel). The ratio of cross sections between the normal (NH) and inverted (IH) neutrino mass 
hierarchies is shown in the lower panel. 
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E [MeV] 


FIG. 4: Same as Fig. [3l but for the i/e-^^Fe reaction. 
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FIG. 5: Same as Fig. [3l but for the t'e-^°®Pb reaction. 
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A2 
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FIG. 6: The ratios race of the number of the detector events induced in mineral oil (CH 2 ), water (H 2 O), and for the 

incoming (anti)neutrino fluxes of the accretion phase in normal (NH) and inverted (IH) mass hierarchy: (a) 


(b) (^) ^e+n*" The three columns correspond to (i) {E°J = 8 MeV, (Flgj = 11 MeV, (E°J = 13 MeV; (ii) 

{E°J = 10 MeV. (El) = 13 MeV, {E°J = 15 MeV; and (iii) {E°J = 12 MeV, {El) = 15 MeV. {E^J = 19 MeV. Al, A2, 
A3 denote the combinations of the luminosity ratios for the accretion phase given in Tab. H] 
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FIG. 7: The ratios rcooi of the number of the detector events induced in mineral oil (CH 2 ), water (H 2 O), and for the 

incoming (anti)neutrino fluxes of the cooling phase in normal (NH) and inverted (IH) mass hierarchy for three configurations 
of initial v average energies. The same notation as in Fig. [6] applies. Cl, C2, C3 denote the combinations of the luminosity 
ratios for the cooling phase given in Tab. |T] 

















